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SUPPLEMENTARY MATERIAL 1 

ON THE EVOLUTION OF MIGRATION IN CHANGING ENVIRONMENTS 

Oana Carja, Uri Liberman, and Marcus W. Feldman 

Consider an infinite, haploid population divided into two denies E x and E y , with 
different selection regimes. Individuals in this population are characterized by two biallelic 
loci, a major locus A/ a and a modifier locus B/b, where the major locus controls the 
phenotype of an individual and the fitness of its phenotype, while the modifier locus is 
assumed to be selectively neutral and controls the migration rate between the two demes. 

We study the evolution of the modifier locus B/b and determine the evolutionarily 
stable migration rate as a function of the pattern of fluctuation in selection experienced 
by the population. We track the frequencies of the four genotypes AB, Ab, aB and ab 
in the population. At each generation, there is recombination and selection in each deme 
separately, after which individuals may migrate between the two demes. As in general 
analyses of neutral modifiers (see Feldman and Liberman 1986), we frame the question in 
terms of the stability of the fixation equilibrium with only B present in the population, 
producing migration rate m^, to invasion by allele b, which produces migration rate mj,. 
We assume these rates to be the same from E x to E y and from E y to E x . This symmetry 
assumption makes our analysis more tractable but relaxing this assumption should not 
change the general conclusions of the model. 

We first present the case when selection is constant in time in both demes; this will 
serve as our reference model. We then consider regimes in which there is fluctuating 
selection within each of the demes. In particular, we ask how the rate of environmental 
volatility affects the stable migration rate. 
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CONSTANT SELECTION 



Let the fitnesses and the frequencies of the four genotypes in each deme be 

Deme E x E y 

Genotype AB Ab aB ab AB Ab aB ab 



Fitness w Xl w X2 w X3 w X4 w yi w V2 w m w yA 



(51.1) 



Frequency x\ x 2 x 3 x 4 yi y 2 yz 2/4- 

Genotypes carrying the B allele at the modifier locus migrate from one deme to the other 
deme at rate rag, while the migration rate is m b when b is present at the modifier locus. 
The recombination rate between the two loci is r in both demes. 

After recombination, selection, and migration, in that order, the change in frequency 
of the genotypes (x[ , x' 2 , x' 3 , x' A and 2/25 2/3? 2/4) i n the two demes from one generation to 
the next is given by 

w Xl (xi - rD x ) w yi (yi - rD y ) 



M x x[ = (1 - m B ) — — + m B 



W X Wy 



M x x> 2 = (1 - m b ) W ^ X2 + TDx) + m h W ^ V2 + rDy) 



A , 1 f, w X3 (x 3 + rD x ) w y3 (y 3 +rDy) 
M x x 3 = (1 - m B ) — Ym B - 



(51.2a) 



W X Wy 



M x x' 4 = (1 - m b ) W ^- rDx) + m b W ^ ~ rDy) 



W X Wy 



and 



A/f 1 ft \ w vi(Vi- rD v) , w Xl {xi-rD x ) 
M y y 1 = (1 - tub) — — + m B - 



w„ w 3 



A/f 1 f, ^ w yAV2 + rD y ) w X2 (x 2 + rD x ) 

M y y 2 = (1 - m b ) — — + m b 



Wy w x 



/ ,t ^ w vs(y3 + rD y ) w X3 (x 3 + rD, 

M y y 3 = (1 - m B ) — — + m B 



(51.26) 



Wy W X 



M y y' 4 = (1 - m b ) W ^- rDy) + m h Wx ^ ~ rD 



W v Wn 



with 



D x = XiX 4 - x 2 x 3 D y = y 1 y 4 -y 2 y 3 (51.3) 

4 4 
w x = ^w Xi Xi w y = ^2wy t yi (51.4) 



w yi ~ w yi 



W U3 ~ W V4- 



(51.5) 



Equation (SI. 5) is due to the assumption that the modifier locus is selectively neutral. 
M x and M y are normalizing factors such that J2i=i x '% = 1 an d Si=i 2/^ = 1- 

We are able to derive closed-form mathematical results if we assume the following 
symmetry relations among the fitness parameters: 

W X1 = W X2 = Wy a = Wy 4 = (l + S) 



Wx 3 = W XA = W yi = W y2 = 1. 



(51.6) 



Thus, if we assume that s > 0, phenotype A is preferred in deme E x , while phenotype 
a is preferred in deme E y and there is symmetry in these fitness differences between the 
demes. 

Boundary equilibria with symmetric selection 

Suppose that initially only B is present at the modifier locus (x2 = X4 = 1/2 = y4 — 0). 
We use the notations x\ = x and ys = z; therefore, X3 = 1 — x, y\ = 1 — z The change in 
the frequencies on the boundary where only B is present is given by 



/ /, Al + s)x 1-z 
x = 1-mB Vm B 

w x w z 

, (l + s)z 1-x 
z =(l-m B ) Ym B , 

w z w x 

with 

w x = (1 + s)x + (1 — x) = 1 + sx 

w y = (1 + s)z + (1 — z) = 1 + sz. 
A simple computation shows that 



(51.7) 



(51.8) 



z'-x' = - ^ Az-x). (51.9) 

(1 + sx)(l + sz) K ' v ; 



Hence, at equilibrium, either z = x or 

1 + a 



= 1. (51.10) 



(l + sx)(l + sz) 
Let 

x = (x 1 ,X2,x 3 ,x 4: ) y = (2/1,2/2,2/3,2/4)- (51.11) 

When 2 = x, at equilibrium, we have the following result, similar to that in Balkau and 
Feldman (1973): 
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Result 1. 

1. On the boundary where only the B allele is present, there is a unique symmetric 
equilibrium (x*,y*) 

x* = (x*,0,l -x*,0) y* = (1 -x*,0,x*,0). (51.12) 

Here, x* is the unique positive root of the quadratic equation 

Q(x) = sx 2 + [(a + 2)m B - s] x + m B = 0. (51.13) 

2. (x*, y*) is internally stable on the boundary where only the B allele is present. 

3. When s > 0 and 0 < mg < 1, then x* > -js. If, in addition, 0 < m B < \ then 



x*>±. 



Proof of Result 1. 



1. When z = x at equilibrium, the mean fitnesses in the two demes are equal, w x = 
1 + sx = 1 + sz = Wy, and the equilibrium equation resulting from (SI. 7) reduces to 
the quadratic equation (SI. 13). Moreover, as 0 < m B < 1, 

Q(0) = -m B <0, Q(l) = (s + l)m B > 0, Q(±oo) > 0. (51.14) 

Therefore there is a unique root x* of (SI. 13) with 0 < x* < 1. 

2. "Near" (x*,y*) and on the boundary where only B is present, (z — x) is small and 
from (SI. 9) its change in magnitude, and thus the internal stability of (x*,y*), are 
determined by the positive factor 

C* = ^ r. (51.15) 

(1 + sx*)(l + sz*) v ' 

(x*,y*) is internally stable if C* < 1, which, since x* = z* and s > 0, is the case if 
and only if s (x*) + 2x* > 1. Using the equilibrium equation (SI. 13), 

s (x*) 2 + 2x* = (1 - m B )(s + 2)x* + m B . (51.16) 
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Thus, as 0 < ins < 1, (SI. 3) is greater than 1 if and only if x* > Since Q(l) > 0 
and Q(x*) = 0, we have x* > if and only if Q^^^) < 0. Indeed, when s > 0, 



Q(T^)=-^y<0, (51.17) 



1_\ _ s(s + l] 
s + 2j~ ( S + 2Y 

as desired. Therefore x* > C* < 1, and (x*,y*) is internally stable. 

3. When s > 0 and 0 < tub < 1, x* > In addition 

OQ) =-|(l-2m B ). (51.18) 
Hence, if 0 < m B < |, Q(|) < 0 and, as Q(l) > 0, we have x* > \. 



On the boundary when only B is present, we may also have equilibria where z / i 
and (SI. 10) is satisfied. These are asymmetric equilibria, and the following result holds: 

Result 2. On the boundary where only the B allele is present, there are two asymmetric 
equilibria (x, y) and (x,y), corresponding to the fixations of either AB or aB, 

x = y = (1,0,0,0), x = y = (0,0,1,0). (51.19) 

Both equilibria are internally unstable. 

Proof of Result 2. If an asymmetric equilibrium with z ^ x exists, then from (SI. 10), 
we have 

1 + sz = 1 + S , (51.20) 

1 + sx v ' 



and therefore 



z = ^\ l-z= S 4±^-. (S1.2V 
s(l + sx) s(l + sx) v 



Substituting these relations into the equilibrium equations resulting from (SI. 7), x satisfies 
the equation 

T(x) = (1 -m)s 2 x(l -x) = 0. (51.22) 

Thus, x = 0 (and z = 1) or x = 1 (and z = 0), giving in the two fixations in either AB 
or aB. We check the internal stability of the fixation in AB, for example. "Near" this 
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fixation, w = 1 — x and z are small, and up to non-linear terms in w and z, their change 
in frequency is 



w = — w + m B [l + s)z 



z = 



1 + 8 

m B 



(51.23) 



l + s 



w + (1 - m B )(l + s)z. 



The characteristic polynomial P(X) associated with the transformation (SI. 23) is 

1 



P(X) = X 2 - (1 - m B ) 



l + s) 



l + s 



A + (1 - m B ), 



(51.24) 



with 



P(l) = 



g 2 (l - rrag) 
l + s 



<0, 



P(+oo) > 0. 



(51.25) 



Thus, -P(A) has a root larger than 1, and this fixation in AB is internally unstable. 
Similarly, the fixation in aB is also internally unstable. 



External stability of the symmetric equilibrium 

Results 1 and 2 entail that, on the boundary where only the modifier allele B 
is present, the unique symmetric equilibrium (x*,y*) is the only internally stable equi- 
librium. It is therefore interesting to explore when (x*,y*) is externally stable to the 
introduction of a modifier allele b determining the migration rate m& in the population. 

The external local stability of the symmetric equilibrium (x*,y*) is determined by 
the linear approximation L*, to the transformation (SI. 2) near (x*,y*). Up to non-linear 
terms, we have 

(?)-(?) + (J)-(?) +t -(s)- 

Here e = (ei, 62, £3, £4) and 5 = (61, 82, S3, 84) are such that £j and <5j are small in magni- 
tude, for i = 1,2,3,4, J2i £ i = J2i^i = °- Here > £ ' = (£1, £2, £3, £4) and 5' = (61,82,63,64) 
have the same properties as £ and 8, respectively. As x\ , y* > 0 for i = 1, 3 and x* = y* = 0 
for i = 2,4, L* has the following block structure 
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ei £ 3 Si S 3 e 2 £4 S 2 



L* = 



0 



V 



Si 
Ss 
£2 
s 4 
S 2 
S4 



(51.27) 



/ 



Here 0 is the 4x4 zero matrix and * is a 4 x 4 matrix that does not affect the eigenvalues 
of L*, which are those of L; n and L ex . Li n determines the internal stability of (x*,y*) on 
the boundary where only the genotypes AB and aB are present. As (x*,y*) is internally 
stable, all eigenvalues of Li n are less than 1 in magnitude. L ex corresponds to the linear 
approximation of the change in genotype frequencies for the genotypes Ab and ab "near" 
(x*,y*). L ex thus determines the external stability of (x*,y*). 

In what follows, we characterize L ex and its largest eigenvalue. The relative relation 
of the largest eigenvalue to 1 gives information about the stability of the equilibrium. 

A straightforward computation shows that 

/(l-m b )A* (l-m b )B* m b C* 



L ex — 



(l-m b )D* (l-m b )C* m b B* 
m b A* m b B* (1 - m b )C* 

m b D* m b C* (l-m b )B* 



m b D* 
m b A* 
(l-m b )D* 
(l-m b )A* ) 



(51.28) 



where 



(51.29) 



(1 + sx*)A* = (1 + s) [1 - r (1 - x*)} 

(l + sx*)B* = (l + s)rx* 

(l + sx*)C* = l-rx* 

(l + sx*)D* =r(l-x*). 
Since 0 < x* < 1 and 0 < r < 1, 0 < m b < 1, L ex is a positive matrix, and by the 
Perron- Frobenius theorem, its largest eigenvalue is positive. Let 5(A) = det (L ex — AI), 
where I is the 4x4 identity matrix, be the characteristic polynomial of L ex . Due to the 
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structure (SI. 28) of L ex (also see Supplementary Material 1 — Appendix A), 5(A) 
factors into the product of two quadratic polynomials Si (A) and S 2 (A): 



S(X) = S 1 (X)-S 2 (X), (51.30) 

where 

Si (A) = A 2 - A [(1 - m h ) (A* + C*) + m b (B* + £>*)] + (1 - 2m b )(A*C* - B*D*) 

S 2 (A) = A 2 - A [(1 - m h ) (A* + C*) - m b (B* + £>*)] + (1 - 2m 6 )(A*C* - 5*ZT). 

(51.31) 

The discriminant of 5i(A) = 0 is 
[(1 - m b ){A* - C*) + m b {B* - D*)] 2 + 4m 6 (l - m 6 )(,4*£>* + 5*C*) 

(51.32) 

+ 4 (1 - m b f B*D* + 4m*,A*C*, 
which is positive since A* , B* , C* , D* are all positive and 0 < m b < 1. In addition, 

(1 + sx*) 2 (A*C* - B*D*) = (1 + s)(l - r). (51.33) 

Thus, when 0 < m b < \ and 0 < r < 1, the two roots of Si (A) = 0 are positive, and they 
are both less than one if and only if Si (A) > 0 and 5^(1) > 0. 

It is not obvious that the roots of S2(A) = 0 are real. However, as the largest 
eigenvalue in magnitude of L* x is positive, in order for it to be less than 1, it is sufficient 
that both positive roots of 52(A) are less than 1. This is the case if S^(l) > 0 and 
5 2 (1) > 0. As 0 < m b < 1 and A*,B*,C*,D* are all positive, we have 

S 2 (l) > Si(l), 5 2 (1) > Si(l), (S1.34) 

so if 5i(l) > 0 and Si(l) > 0, it is also true that S 2 (l) > 0 and S 2 (l) > 0. As (x*,y*) is 
internally stable, from Result 1 we know that (1 + sx*) 2 > (1 + s), and so, by (SI. 33), 
the constant term of both Si (A) and S 2 (A) satisfies 

0 < (1 -2m b )(A*C* -B*D*) < 1, (51.35) 

provided 0 < r < 1 and 0 < m b < \. It is therefore impossible that the two positive 
roots of Si (A) = 0 are both larger than 1. Hence, they are both less than 1 if and only if 
5i(l) > 0, in which case also S^(l) > 0. 
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To sum up, (x*,y*) is externally stable if and only if 5i(l) > 0. Computing 5i(l), 
we find that 5i(l) > 0 if and only if 

(m b -m B (l + sx*y 2 (1 - r)s [x*(s + 2) - 1] > 0. (51.36) 

As x* > 0 < r < 1 and s > 0, (SI. 36) holds if and only if mj, > mg. We have proved 

Result 3. The internally unique stable symmetric equilibrium is externally stable towards 
the introduction of a modifier allele b provided its associated migration rate nib is larger 
than mB- 

This result is in complete accordance with the reduction principle of Feldman and 
Liberman (1986) for modifiers of migration rates. It implies that the uninvadable stable 
migration rate is zero. We now study the case when selection pressures are no longer 
constant, but fluctuate in time. 



PERIODICALLY FLUCTUATING SELECTION 



Suppose now that, within each deme, the selection regime varies temporally. We assume 
two possible types of selection regimes. Under T\ and T<i- In Ti, the fitness coefficients 
are assumed to be 



while in they are 



Deme 
Allele 

Fitness 

Deme 
Allele 
Fitness 



E x 
A a 
1 + s 1 



Ey 

A a 
1 1 + s, 



(51.37) 



E x E. 
A a 
1 1 + s 



y 

A a 
+ s 1. 



(51.38) 



For simplicity, we assume that within each deme, phenotype A is favored in one 
temporal state and phenotype a is favored in the other. The environments in which allele 
A does better are T\ in deme E x and T-i in deme E y . The selection regimes fluctuate 
between T\ and T 2 . We first assume periodically fluctuating selection and explore the 
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evolution of the migration rate, as a function of the length of the period of the fluctuating 
selection, which we also refer to as the environmental period. In the general case, an 
environmental cycle consists of t\ + t<i selection steps, the first t\ of type T\ followed 
by T2 of type T 2 selection. We are able to derive closed-form solutions and fully explore 
the dynamics of the system in the case t\ = T2 = 1, i.e. the environment changes 
every generation. For analytical tractability, we further need to assume that all selection 
coefficients are equal. 

Let T be transformation of population state given in (SI. 2) with associated fitness 
parameters of Ti, and let T be the transformation associated with T2 selection. The 
transformation of the population state after the (t\ + T2) phases of selection is given by 

= rofo---oT o foro---or f*^ . (51.39) 

T2-times Ti-times 

We start the analysis with the case t\ = T2 = 1, namely 

x' \ ~ = / X 



y'J =ToT U (5L40) 



which can be written as a two-phase transformation, 



✓ )- f (P)' W- T {y)- (SL41) 



Boundary equilibria with symmetric selection 

Assume that, initially, only B is present at the modifier locus, and let x be the 
frequency of AB in deme E x and z the frequency of aB in deme E y . Then x = (x, 0, 1 — 
x, 0), y = (1 — z, 0, z, 0), and using (SI. 2), we obtain 



and 



~ Al + s)x 1 
x = (1 - m B )— h m B - 



1 + sx 1 + sz 

,i \ (1 + s ) z 1 — X 

z = (1 - m B )— h 



1 + sx ' 



X =(1 - mB) l + S (l-5) +mB l + 

, 5 (l + s)(l-x) 

^ =(1 - mg) l + g (l-^) +m " 1 + 8 (1-X) 
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(51.42) 



(51.43) 



Then, 

(~ ~\ 1 + 5 ( \ 

(Z — X) = — rr- r{Z — X) 

K ' (1 + sx)(l + sz) y ' 



(51.44) 



Hence, 



<z'-x') = - . r + ^ — -, rrU-x). (51.45) 

V ; (l + sx)(l + sz)[l + s(l-3:)}[l + s(l-z)] K ' V ; 



Thus, at equilibrium, we either have z = x or 



(1 + S)2 =1. (51.46) 



(1 + 8X)(1 + SZ) [1 + S(l - X)} [1 + S(l - 5)] 

When z = i, we have a symmetric equilibrium characterized by the following result: 
Result 4. 

1. On the boundary where only the modifier allele B is present, if 0 < ms < 1 and 
s > 0, a unique symmetric equilibrium (x, y) exists with 

x = (x,0, 1 -x,0), y = (1 -x,0,x,0). (51.47) 

Here, x is the unique positive root of R(x) = 0, where 

R(x) = sx 2 + [2 - m B (s + 2)] x - (1 - m B ). (51.48) 

2. When 0 < ms < 1, x > If, in addition, 0 < vtlb < \, then 0 < x < \. 

3. (x, y) is internally stable on this boundary 

Proof of Result 4. 

1. At a symmetric equilibrium, z = x and therefore z = x. From (SI. 42) and (SI. 43) 

~ = [(1 - m B )(s + 1) - m B ] x + m B 

1 + sx { ' 

_ [(1 - tub) - m B {s + 1)] x + m B (l + s) 
X_ l + s(l-x) ' ^ ' 

Substituting (S1.49) into (S1.50) gives 

(s + 2)m B {sx 2 + [2 - m B (s + 2)] x - (1 - m B )} = 0. (51.51) 
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As rriB > 0 and s > 0, x satisfies the quadratic equation R(0) = 0. As i?(±oo) > 0, 
R(x) = 0 has two real roots, one positive and one negative. But, as 0 < mg < 1, 



R(l) = (s+ 1)(1 -m B ) > 0, (51.52) 



and thus R(x) = 0 has a unique root x with 0 < x < 1. 
2. Observe that 



As .R(O) < 0 when 0 < m B < 1 and R(l) > 0, when 0 < m B < 1, we have x > 
If, in addition, 0 < tub < \, we have 0 < x < ^. 

3. In view of (SI. 45), the symmetric equilibrium (x, y) is internally stable provided 

o < !■ (^1-54) 

(l + sx) 2 [l + s(l-5)] 2 

Using (S1.49), 

(1 + sx)[l + s(l - 5)] = (1 + s) + sm B [(s + 2)x - 1]. (51.55) 

When 0 < mg < 1, x > an d therefore (1 + sx)[l + s(l — x)] > 1 + s, and thus 
(SI. 54) is satisfied. 

As in the constant environment case, it is possible that asymmetric equilibria exist 
with z 7^ x and (SI. 46) is satisfied. For example, the two fixations in AB or in aB, in 
both demes, are such equilibria. They correspond to i = 1, z = 0 or i = 0, z = 1, 
respectively. We can show that both fixations are internally unstable. Our computer 
simulations suggest that there are no stable asymmetric equilibria. 



External stability of the symmetric equilibrium 

We next check the external stability of the internally stable symmetric equilibrium 
(x,y) to the introduction of b at the modifier locus, changing the migration rate from mg 
to mi,. As the transformation of the population state is T o T, following the same type 
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of analysis as in the constant environment case, the matrix L ex determining the 
stability of (x, y) is given by 



L ex — L ex o L 



'ex? 



where, as in (SI. 28), 

/(l-m b )A (l-m b )B 



L ex — 



m b D \ 
rribA 



m b C 

1 — m b )C m b B 
m b B (1 — m b )C (1 — m b )D 
m b C (1 — m b )B 



J ex 



m b A 
\ m b D 
As in (SI. 29) we have 



l-m b )D 
m b A 
\ m b D 

/(l-m b )A (l-m b )B m b C 

(1 — m b )D (1 — m b )C m b B 

m b B (1 — m b )C (1 — m b )D 

m b C (l-m b )B (l-m b )Aj 



(l-m b )A/ 

m b D \ 

m b A 



(1 + sx)A = (1 + s) [1 - r (1 - x)] 



■ sx) 



)B = (1 + s)rx 
(1 + sx)C = l- rx 
(1 + sx)D = r(l - x), 



x 



(l + sx)A = (1 + s) [1 -r(l 
(1 + sx)S = (1 + s)rx 
(1 + sx)C = l- rx 
(1 + sx)D = r(l - x). 
The relation between the equilibrium point (x, y) and (x, y) is such that 



(x,y) ^> (x,y) ^> (x,y). 



We can write 



(a e h d\ 

b f 9 c 

c 9 f b 

\d h e a J 
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(51.63) 



where 

a = (1 - m b ) 2 AA + (1 - m b ) 2 BD + m\CA + m 2 b DD 

b = (l- m h f DA + (1 — m h f CD + m 2 b BA + ml AD 

c = m b (l-m b ) AA + BD + CA + DD 

d = m b (l-m b ) DA + CD + BA + AD 

e = (1 - m b f AB + (1- m b f BC + m\CB + m 2 b DC 

/ = (1 - m 6 ) 2 55 + (1 - m 6 ) 2 CC + m^B + m^lc 

g = m b (l-m b ) AB + BC + CB + DC 

h = m b (l - m b ) \b~B + CC + BB + AC . 

From (SI. 57) and (SI. 58), the characteristic polynomial -D(A) = det (L ex — AI) of L ex , 
factors into the product of two quadratic polynomials: -D(A) = D 1 (\)D2(\). In fact, 

/ a + d- A e + h 0 0 \ 

b+c f+g-X 0 0 

e g f - g- A 6-c 

d /i e — /i a — d — X/ 



D(X) = 



\ 



(51.64) 



with 



D 1 (X) = X 2 -(a + d + f + g)X + (a + d)(f + g) - (b + c)(e + /i) 
D 2 (A) = A 2 - (a - d + / - g)X + (a - d)(f - g) - (b - c)(e - h). 



(51.65) 



Based on Supplementary 1 — Appendix A, the constant terms of -Di(A) and -D2(A) 
are both equal to 

(1 - 2m b f (1 + sx)~ 2 [l + s(l- x)] 2 (1 + sf (1 - rf . (51.66) 

These constant terms are positive for 0 < r < 1 and 0 < m b < | and are less than 1 due 
to the internal stability of (x,y). Observe that the discriminant of -Di(A) = 0 is 



(a + d + f + g) 2 - 4(la + d)(f + g)-(b+c)(e + h)} 



(567) 



which reduces to 



[(a + d)-(f + g)] 2 + 4(6 + c)(e + h)> 0, 
15 



(568) 



since a, b, c, d, e, /, g, h are all positive. Therefore -Di(A) = 0 has two positive roots whose 
product is less than 1. These positive roots are less than 1 provided both D 1 (l) and D[(l) 
are positive. As it is impossible that both roots are larger than 1 (their product is less 
than 1), we only need the condition -Di(l) > 0. Since the matrix L ex is positive, we know 
that its largest eigenvalue is positive. For the external stability of (x, y), we also need to 
prove that when the roots of D 2 (\) = 0 are real and positive (their product is positive) 
they are both less than 1. This is the case when -D2(l) > 0 and D' 2 (l) > 0. It is clear that 



D 2 (1)>D 1 (1), D' 2 (l) > D[(l). 



(51.69) 



Therefore, when D x (l) > 0 and D[(l) > 0, this also implies D 2 (l) > 0 and D' 2 (l) > 0. 
To sum up, (x, y) is externally stable provided -Di(l) > 0 where 



D 1 (l) = 1 - (a + d + / + g) + (a + d)(f + g) - (b + c)(e + h), 

and (a + d)(f + g) - (b + c)(e + h) is given in (SI. 66). 
Let us compute (a + d + f + g). Following (SI. 63), 

(a + d + f + g) = \(1 - m b f A + mlC + m b (l - m b )(B + D) 



+ 



+ 



(1 - m b f C + m 2 b A + m b (l -m b )(B + D) 
(1 - m b f B + m 2 b D + m b (l - m b ) {A + C) 
(1 - m b f D + m 2 b B + m b (l - m b ) (A + C) 



A 
C 
D 
B 



which can be rewritten as 



(a + d + f + g) =(l-m b ) 



AA + BD + CC + DB 
+ m b (l - m b ) \(B + D) (A + C) + (A + C) (B + D) 
AC + BB + CA + DD 



+ m\ 



Thus, 



(a + d + f + g) = (l-2m b ) 
+ m b 



AA + BD + CC + DB 
(B + D)(A + C) + (A + C)(B + D) 
+ m 2 b (A + C - B - D) (A + C + B -D). 
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(51.70) 



(51.71) 



(51.72) 



(51.73) 



By (S1.59) and (S1.60), 

(1 + sx) [1 + s(l - x)] [AA + BD + CC + DB] 



(51.74) 

= (1 + s) [2(1 - r) + r 2 ] + r 2 s(s + l)x - r 2 sx - r 2 s 2 xx, 



(1 + sx) [1+8(1 - x)} [(B + D) (A + C) + (A + C) (B + C) 

= 2r 2 (l + sx) [1 + s(l - x)] + r(l - r)(s + 2) [(s + 2) + s(x + x)] , 
(1 + sx) [1 + s(l - x)] [A + C - B - D] [A + C - B - D] = 

= (s + 2) 2 (l-r) 2 . 

Using (SI. 66), Di(l) is equal to 

(1 + s) [2(1 — r) + r 2 l + r 2 s(s + l)x — r 2 sx, — r 2 s 2 xx 
1 " (1 " 2mb) (! + .*)[! + .(!-*)] 

2r 2 (l + sx)[l + s(l - x)] + r(l - r)(s + 2)[(s + 2) + s(x - x)] 



(51.75) 



(51.76) 



- m b 



+ 



(1 + sx)[l + s(l - x)} 

2 (s + 2) 2 (l-r) 2 
(1 + sx)[l + s(l - x)} 

(1 - 2m b f (s + if (1 - rf 



(51.77) 



m 



(1 + sx) 2 [1 + s (1 - x)} 2 



We next show that -Di(l) can be represented as 

D 1 (l) = (l-r)s(m B -m b )A(r), (51.78) 

where A(r) is a linear function of r that is positive for all 0 < r < 1. Therefore the sign 
of -Di(l) coincides with that of (ms — mb), and -Di(l) > 0 when mg > m b . We thus have 
the following result: 

Result 5. The unique internally stable symmetric equilibrium (x, y) is externally stable 
to the introduction of the allele b at the modifier locus provided vtib > mb- Thus, higher 
migration rates are favored and the stable uninvadable migration rate is 1. 

Computation of -Di(l). 
Observe that 

r 2 [(l + s) + s(s + l)x -sx- s 2 xx] = r 2 (l + sx)[l + s(l - x)], (51.79) 
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so (SI. 77) simplifies to 

2 2(l + a)(l-r) r (i- r )( g + 2 )[( g + 2) + g (x-x)] 

1 b) (1 + sx)[l + 8(1 - x)] mb (l + sx)[l + s(l-x)} 



(s + 2) 2 (1 - rf (1 - 2m b f (s + l) 2 (1 - rf 



" >,b (1 + 8X)[1 + S(l - x)] ' (1 + S x) 2 [1 + 8 (1 - x)f _ 



(51.80) 



Clearly -Di(l) of (SI. 80) has a factor of (1 — r), and in fact 

Di(l) = (l-r)/(r), (51.81) 

where f(r) is a linear function of r, for 0 < r < 1. Now 

, m o o 2(l + g) ( s + 2 )[( s + 2) + s(x-x)] 

/(l) = 2-(l-2m 6 ) (1 + a _ )[1 + a(1 _ x)] -m b (1 - - ^ _ „ }] . (51.82) 

Following (SI. 55) we have 

(1 + sx)[l + s(l - 5)] = (1 + s) + sm B [(s + 2)x - 1]. (51.83) 

We also have an equivalent expression for (SI. 83) in terms of x, namely 

(1 + sx)[l + s(l - 5)] = (1 + s) + sm B [(s + 1) - (s + 2)5]. (51.84) 

Also, whereas x > we have x < Applying all of this to (SI. 82) and using the 

fact that 

(1 + sx)[l + s(l - 5)] = (1 + s) + + (s + 2)(x - x)], (51.85) 

we get that 

(1 + sx)[l + s(l - x)]/(l) = 2(s + 1) + sm B [s + (s + 2)(x - 5)] 

- 2(1 - 2m b )(s + 1) - m 6 (s + 2)[(s + 2) + s(x - x)] 
= sm B — mb (s + 2) — 4(s + 1) + s(s + 2)(m B — m^ix — x) 
= s 2 {m B - m b ) + s(s + 2)(m B - m b )(x - x). 



(51.86) 

Thus 

(1 + sx)[\ + s(l - x)]/(l) = s{m B - m h )[s + (s + 2)(x - x)]. (51.87) 
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But as (s + 2)x > 1, (a + 2)x < (s + 1), 

s + (s + 2)(x-x) = [(a + 2)x - 1] + [(a + 1) - (a + 2)x] > 0. (51.88) 

It follows that the sign of /(l) is the same as the sign of (mg — m^). 
We now compute /(0): 



2(1 + a) 2 (a + 2) 



2 



/(0) = 1 - (1 - 2m 6 ) - ■ ■ ' - m 2 



+ 



(1 + Sx)[l + 8(1 ~ X)} b (1 + SX)[1 + S(l - X)} 

(l-2m 6 ) 2 (s + l) 2 



(51.89) 



(l + sx) 2 [l + s(l-x)] 2 ' 



Using the expression (SI. 83) for the product of the two mean fitnesses, we get 
(1 + sxf [1 + s (1 - x)} 2 /(0) = {(1 + s) + sm B [(s + 2)x - l]} 2 

- 2(1 - 2m b )(s + 1) {(1 + s) + sm B [(a + 2)x - 1]} 

- m^(s + 2) {(1 + s) + sm B [(a + 2)x - 1]} 
+ (l-2m 6 ) 2 (s + l) 2 . 

(51.90) 

In (SI. 90) we replace the x 2 term using the equilibrium equation (SI. 48) to give 

(1 + sxf [1 + a(l - x)f f(0) = (m B - m b )s{m b s(s + 1) 

- m% (a + 2) 2 + m B Us + l)(a + 4) - m b (a + 2) 2 

+ m B (s + 2)x m B (s + 2) 2 + m 6 (s + 2) 2 -4(s + l) }. 

(51.91) 

The right-hand side of (SI. 91) is (mg — m b )s multiplied by 

m b s{s + l)+m B (s + 2) 2 (m B + m b )[x(s + 2) - 1] + m B (s+ l)[(s + 4) -4x(a + 2)]. (51.92) 

We will show that (SI. 92) is always positive. In fact, (SI. 92) is equal to 

m b s(s + 1) + m B ■ m b (a + 2) 2 [x(s + 2) - ll + m% (a + 2) 2 [x(s + 2) - ll 

(51.93) 

+ m B (s + 1) [(a + 4) - Ax(s + 2)] . 
From the equilibrium equation (SI. 48) we get that 

m B [(s + 2)x - 1] = sx 2 + 2x-l. (51.94) 
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Hence (SI. 93) is equal to 



m b s(s + l) + m B -m b (s + 2) 2 [x(s + 2) - l] + 

+ m B (a + 2) 2 [sx 2 + 2x - l] + m B {s + 1) [(s + 4) - 4x(s + 2)] . 



(51.95) 



The last two terms have a factor mg that multiplies 



(a + l)(s + 4) - (a + 2) 2 + (s + 2) 2 x(2 + sx) - 4x(s + l)(s + 2) 
= s + (s + 2)x [(a + 2)(2 + sx) - 4(a + 1)] 
= s + (s + 2)x [(a + 2)sx - 2s] 

= a [(a + 2) 2 x 2 - 2(s + 2)x + lj = s [(s + 2)x - l] 2 , 



(51.96) 



which is positive. To sum up, /(0) also has the same sign of (mg — mj), and so 



Di(l) = (l-r)s(m B -m b )A(r), 



(51.97) 



where A(r) is a linear function of r that is positive for all 0 < r < 1. 

We saw that in a temporally constant environment the reduction principle holds and 
smaller migration rates are always favored. On the other hand, with fluctuating selection of 
period 2, non-zero migration rates evolve and the reverse of the reduction principle holds: 
a higher migration rate can always invade in a population fixed on a smaller resident 
one and the evolutionarily stable migration rate is 1. An important result is that this 
evolutionarily stable switching rate maximizes mean fitness at equilibrium. This result 
is deeply rooted in classic population genetics theory and is connected to early work on 
neutral modifiers and the mean fitness principle of Karlin and McGregor (1972, 1974). 

The mean fitness and external stability. 

We saw that in a constant environment the symmetric equilibrium is externally stable 
if nib > itib and unstable if m b < m B , so that smaller migration rates are favored. On 
the other hand, with fluctuating selection of period 2, non-zero migration rates evolve; the 
symmetric equilibrium is stable if m B > m& and unstable if m B < mt,. This phenomenon 
is related to the behavior of the mean fitness at equilibrium in the following way. 
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Result 6. The mean fitness at the symmetric equilibrium is 

i. a decreasing function of tub in a constant environment, 

ii. an increasing function ofrriB in a period 2 cycling environment. 

Proof of Result 6. 

i. In a constant environment, the mean fitness w* at the symmetric equilibrium (x*, y*) 
is w* = 1 + sx*, and it is a decreasing function of if is negative. Using the 
equilibrium equation, 

dx* 1 - (s + 2)x* 



dm B 2sx* + [(s + 2)m B - s] ' 
As x* > in order for J^-^ to be negative, it is sufficient that 



(51.98) 



x* > s ~ m f + 2 \ (51.99) 

This follows easily from the fact that Q(x) in (SI. 13) satisfies Q(0) < 0, Q(x*) = 0, 
and Q ( s ' m f s (s+2) ) < 0. 

ii. With a fitness cycle of period 2, the mean fitness w at the symmetric equilibrium 
(x, y) is 

w = (1 + s) + sm B [(s + 2)x - 1]. (51.100) 
w is an increasing function of if > 0. Now 

-^ = sUs + 2)x-l]+s(s + 2)m B ^-. (51.101) 
om B om B 

Thus > 0 provided > 0. Using the equilibrium equation R(x) = 0 for x, 

we have 

gg = x(s + l)-l 
<9m B 2sx+ [2-m B (s + 2)] 

Since [x(s + 1) - l] > 0, we conclude that - > 0 if 

_ > ma( . + 2)-2 | (S1 1Q3) 

which follows from i?(0) < 0, i?(x) = 0, and R ( mB(s + 2) ~ 2 ^) < 0. 
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SUPPLEMENTARY MATERIAL 1 - APPENDIX A 



The family of the 4x4 matrices of the form 

mC mD \ 

mB mA 

(l-m)C (l-m)D 

(l-m)B (l-m)A/ 



L = 

V mD 
has the following properties: 
1. If P is the 4x4 matrix 



/(l-m)A (l-m)B 
[l-m)D (l-m)C 
mA mB 



mC 



(1 0 0 1\ 

0 110 

0 0 10 

\0 0 0 1/ 



then 



PLP 



-l 



M 0 

* N 



M 



N = 



where 0 is the 2x2 zero matrix and 

' (l-m)A + mD (l-m)B + mC\ 
mA+ (l-m)D mB + (1 - m)C ) 

(1 — m)C — mB (1 — m)D — mA \ 
(1 — m)B — mC (1 — m)A — mD J 

Hence L is similar, via P, to the block matrix (SI. A3). 
2. The determinants of M and N are equal: 

det(M) = det(N) = (1 - 2m)(AC - BD). 



(S1.A1) 



(S1.A2) 



(S1.A3) 



(S1.A4) 



(S1.A5) 



3. The characteristic polynomial of L factors into the product of the characteristic poly- 
nomials of M and N, namely 



det (L - I 4 ) = det(M - AI 2 ) • det(N - AI 2 ), 

where L is the identity matrix of order % for i = 2, 4. 
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(S1.A6) 



4. If Li,Li2, . . . , L n are matrices of the form (Sl.Al) with corresponding matrices M; 
and Nj as in (SI. A3), for i = 1, 2, . . . , n, then 

n'iiM, o \ 



p p p-' = 



vi=l 



nr=i n, 



(S1.A7) 



Therefore IliLi * s similar, via P, to the block matrix (S1.A7). 

5. The characteristic polynomial of nT=i factors into the characteristic polynomials 

of nr=iM, and nr=iN, 

6. The matrix Yl7=i also has the form 

/a e h d\ 

b f 9 c 

c £ / 6 

V d /i e a / 

and the corresponding matrices Y\7=i anc ^ IliLi have the form 



(S1.A8) 



a + <i e + h 
b+c f+g 



and 



/ -9 b-c 
e — h a — d 



(S1.A9) 



respectively. 

7. If D(\), Di(X), and A (A) are the characteristic polynomials of fliLi ^i, niLi 
and nr=i ^i, respectively, then 

D(A) = A(A)-A(A), 



where 

A (A) = A 2 - (a + d + / + 5) A + (a + d)(/ + g) - (b + c) (e + fc) 

A(A) = A 2 - (a - d + / - <?)A + (a - d)(/ - <?) - (b - c)(e - fc). 
Moreover, the two constant terms of A (A) and A (A) are both equal 

(a + d)(/ + g) - (b + c)(e + h) = (a - d)(f - g) - (b - c)(e - h), 

and (S1.A12) is equal to 

n n n 

11 det(Mi) = J] det(N0 = (1 - 2m)" JJ^Ci - Si A), 

i=l i=l i=l 

where Aj, A,Ci, A, are the parameters associated with Lj, for z = 1,2, 
(Sl.Al). 



(S1.A10) 



(S1.A11) 



(S1.A12) 



(S1.A13) 
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SUPPLEMENTARY MATERIAL 2 



ON THE EVOLUTION OF RECOMBINATION IN CHANGING ENVIRONMENTS 

Oana Carja, Uri Liberman, and Marcus W. Feldman 

Consider an infinite sexual haploid population in which each individual is defined 
by three biallelic loci. The two major loci A/ a and B/b control the phenotype and thus 
the fitness of the individual. The third locus M/m is a modifier locus that controls the 
recombination rate between the two major loci. This modifier locus is otherwise selectively 
neutral. 

We study the evolution of the modifier locus M/m and determine the evolutionarily 
stable recombination rate as a function of the pattern of fluctuation in selection experi- 
enced by the population. As in general analyses of neutral modifiers (see Feldman and 
Liberman 1986), we frame the question in terms of the stability of the fixation equilibrium 
with only M present in the population, to invasion by an allele m, introduced in the pop- 
ulation by a rare mutation. To that end, we track the frequencies of the eight genotypes 
MAB, MAb, MaB, Mab, mAB, mAb, maB, mab. At each generation, the population 
experiences random mating, recombination, and selection. 

As the modifier locus is assumed to be selectively neutral, there are four fitness 
parameters determined by the two major loci 

major loci genotype AB Ab aB ab 

J (S2.1) 
fitness Wi W2 wz W4. 

There are three possible recombination rates depending on the mating type at the 
modifier locus 

modifier locus genotype MM Mm mm 

(52.2) 

recombination rate r 1 r 2 r 3 . 
Suppose that the three loci are ordered with the modifier locus located on one side of 
the two major loci. Let R be the recombination rate between the modifier locus and the 
two major loci, and we assume no interference between recombination events occurring in 
the two intervals separating the two major loci and between the modifier locus and the 
nearest major locus. 
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Let x = (#i, X2, ■ ■ ■ , x 8 ) be the frequency vector of the eight genotypes in the present 
generation, and let x' = (x^x^, . . . ,x' g ) be the frequency vector of the next generation. 
The change in genotype frequencies from one generation to the next is given by 



x\ + R(x2X$ + X3X5 + X4X5 — X1X7 — xixs — xixe) 
+ ri(x 2 x 3 - X1X4) 
+ r 2 (x 2 x 5 + x 2 x 7 - xix 8 - xiXq) 
+ i?r 2 (2xix 6 - 2x 2 x 5 + x 3 x 6 + xix 8 - x 4 x 5 - x 2 x 7 ) 



wx' 2 = W2 



x 2 + R(xixe + x 3 x 6 + x 4 x 6 - x 2 x 5 - x 2 x 7 - x 2 x 8 ) 
+ ri(xix 4 - x 2 x 3 ) 
+ r 2 (xix 6 + xix 8 - x 2 x 5 - x 2 x 7 ) 
+ i?r 2 (2x 2 x 5 - 2xix 6 + x 2 x 7 + x 4 x 5 - xix 8 - x 3 x 6 ) 



wx' 3 = W3 



x 3 + R(xix 7 + x 2 x 7 + x 4 x 7 - x 3 x 6 - x 3 x 8 - x 3 x 5 ) 
+ ri(xix 4 - x 2 x 3 ) 
+ r 2 (x 4 x 7 + x 4 x 5 - x 3 x 8 - x 3 x 6 ) 
+ i?r 2 (xix 8 + x 3 x 6 + 2x 3 x 8 + x 4 x 5 - x 2 x 7 - 2x 4 x 7 ) 



wx' A = W4 



X 4 + R(X\X% + X 2 X 8 + X 3 X 8 — X 4 X5 — X 4 X6 — X 4 X 7 ) 

+ ri(x 2 x 3 - xix 4 ) 
+ r 2 (x3X 6 + x 3 x 8 - x 4 x 5 - x 4 x 7 ) 
+ i?r 2 (x 2 x 7 + x 4 x 5 — 2x 3 x 8 — xix 8 — x 3 X6 + 2x 4 x 7 ) 



wx' 5 = Wi 



x 5 + i?(xix 6 + XiX 7 + XiX 8 - x 2 x 5 - x 3 x 5 - X 4 X 5 ) 
+ r 2 (xix 6 + x 3 x 6 - x 2 x 5 - x 4 x 5 ) 
+ r 3 (x 6 x 7 - x 5 x 8 ) 

+ i?r 2 (2x 2 x 5 - 2xix 6 + x 2 x 7 + x 4 x 5 - Xix 8 - x 3 x 6 ) 
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WXq = w 2 x 6 + R(x 2 x 5 + x 2 x 7 + x 2 x 8 - x 1 x 6 - x 3 x 6 - X 4 X 6 ) 
+ r 2 (x 2 x 5 + x 4 x 5 - xix 6 - x 3 x 6 ) 
+ r 3 (x 5 x 8 - 2:62:7) 

+ Rr 2 (2x 1 x 6 - 2x 2 x$ + Xix 8 + x 3 x 6 - x 2 x 7 - x 4 x 5 ) 

Wx' 7 = W3 X7 + i?(x3X5 + X3X6 + X3X8 — X1X7 — X2X7 — X4X7) 

+ r 2 (xix 8 + x 3 x 8 - x 2 x 7 - X4X7) 
+ ^3(3:52:8 - a: 6 x 7 ) 

+ Rr 2 (2x4x7 - 2x 3 x 8 + X4X5 + x 2 x 7 - xix 8 - x 3 x 6 ) 

WXg = W4 X 8 + i?(x4X5 + X4X6 + X4X7 — XiX 8 — X 2 X 8 — X3X 8 ) 

+ r 2 (x 2 x 7 + X4X7 - xix 8 - x 3 x 8 ) 
+ r 3 (x 6 x 7 - x 5 x 8 ) 

+ i?r 2 (xix 8 + x 3 x 6 - x 2 x 7 - X4X5 + 2x 3 x 8 - 2x4X7) 
Here w is a normalizing factor such that Yli=i x 'i = ^- 



(52.3) 



We first present the case when selection is constant in time in both demes; this will 
serve as our reference model. We then consider regimes in which there is fluctuating 
selection through time. In particular, we ask how the rate of environmental volatility 
affects the stable recombination rate. 
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CONSTANT SELECTION 



We are able to derive closed-form mathematical results if we assume the following 
symmetry relations among the fitness parameters: 

wi = W4, 1U2 = W3. (52.4) 

Also, without loss of generality, we will assume that W\ > W2; that is 

wi = W4 = 1 + s, w 2 = W3 = 1 (52.5) 

with s > 0. 

Equilibria with symmetric selection 

Suppose that initially only allele M is present at the modifier locus. In this case the 
transformation (S2.3) reduces to 

wx'i = wi(xi — r\D) 
wx' 2 = w 2 (x2 + riD) 

(52.6) 

wx' 3 = w 3 (xs + riD) 
wx' 4 = 10^x4 — riD), 
as Xi = x\ = 0 for i = 5, 6, 7, 8. Here 

D = X1X4 — X2X3 

(52.7) 



n 

W 

i=l 



= ^2 WiXi ~ r iD(wi - w 2 - w 3 + W4). 



We now look for equilibria of (S2.6) and check their internal stability. We have the 
following result: 

Result 1. On the boundary where only the M modifier allele is present, the possible 
equilibria are the four fixations and a symmetric equilibrium x* = (x*, | — x*, | — x*, x*) . 
The fixations in AB and in ab are internally stable, whereas the other equilibria are 
unstable. 

Proof of Result 1. 
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As Wi = and W2 = W3, from (S2.6), at equilibrium 



w(xi - X 4 ) = Wl(xi - X 4 ), iu(x 2 - X 3 ) = W 2 (X2 - X 3 ). 



(52.8) 



Hence, x\ = X4 or w = w±, and X2 = X3 or w = 102- If it; = W\, then, from (S2.6) again, 
X\ = Xi — r\D, and as 0 < r\ < 1, D = 0. Similarly, if w = w 2 , then also D = 0. If 
11 7^ 14 and w = Wi, then X2 = X3 as otherwise w = W2, which is impossible. As in this 
case D = 0, we have 

it; = it;i = wi(xi + £4) + u>2(x2 + £3), (52.9) 

and as u>i > u>2, (S2.9) implies that X2 = X3 = 0 and X1+X4 = 1. Since D = X1X4— x 2 X3 = 
0, we then have x\ = 1, £4 = 0 or x\ = 0, £4 = 1. Therefore, the possible equilibria when 
either x\ 7^ x 4 or X2 7^ X3 are the four fixations (1, 0, 0, 0), (0, 1,0, 0), (0, 0, 1,0), (0, 0, 0, 1) 
in AB, Ab, aB, ab, respectively. In addition, if x\ = X4,X2 = £3, there can be an 
equilibrium (x, \ — x, \ — x, x) which satisfies 

wx = Wi(x — r\D) 

, ,1 \ 2 1 

- I - - X = x - - 
2 / 4 



D = x z 



(52.10) 



W = 2(lUi - W2)(l - ri)x + ylfi + \1 - y^2j • 



Thus 



X = 



Wix(l — r\ 



2(wi - w 2 )(l - ri)x + ^-wi + (l - ^) w 2 ' 



(52.11) 



Therefore Q(x) = 0 where 
Q(x) =2(w;i-w; 2 )(l-ri)x 2 



W2 - wi{\ - n] 



ri 



x-^wi. (52.12) 



Now Q(0) = -jWi < 0 and Q (|) = ^u>2 > 0. Hence, Q(x) = 0 has a unique root x* 
satisfying 0 < x* < \ determining an equilibrium x* = (x*, \ — x*, \ — x*,x*). 



We next check the internal stability of these five possible equilibria. The fixation in 
AB is internally stable, as its local stability matrix is 



I l+s 

0 

V 0 



0 
1 

l+s 

0 



ri \ 
l+s \ 

ri 
l+s 

l-rj 



(52.13) 
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with all of its eigenvalues positive and less than 1. Similarly the fixation of ab is stable. 
The local stability matrix associated with the fixation of Ab is 



(1 + 8 nil + s) 
o l - n 
V o n(i + s) 



o \ 

0 

1 + 8/ 



(52.14) 



which has two eigenvalues larger than 1. Similarly, the fixation in aB is also internally 
unstable. 



We next check the stability of x*. At x* we have x\ = £4, X2 = £3. Then "near" x*, 
up to non-linear terms, 

x[-x' 4 =^(x 1 -x 4 ), x' 2 - x' 3 = ^(x 2 - x 3 ). (52.15) 
w w 

Hence x* is internally stable provided w* > wi, w 2 - As w± > W2 and 

w* = 2( Wl - w 2 )(l - n)x* + T -w x + (l - ^ w 2 , (52.16) 

w* > W2- w* > wi if and only if 

2( Wl - w 2 )(l - r)x* > (l - r -) ( Wl - w 2 ), (52.17) 

or, as (wi — ^2) > 0, if and only if 

2(l-riK>l-y, (52.18) 

which is impossible when 0 < r\ < 1 and 0 < x* < \. Therefore x* is not internally 
stable. 

When the environment is constant, the only possible stable equilibria withM fixed are 
the two fixations in AB and ab. These fixations do not depend on n, the recombination 
rate determined by M We now study the case when selection pressures are no longer 
constant, but fluctuate in time. 
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PERIODICALLY FLUCTUATING SELECTION 

Suppose now that the environment varies temporally. We assume two possible types of 
selection regimes such that the fitness parameters change as follows: 



phenotype AB Ab aB ab 

environment 1 wi W2 W3 W4 (52.19) 

environment 2 w\ w>2 W3 W4, 

where we can write 

Wi = w 4 = 1 + s W2 = W3 = 1 

(52.20) 

Wi = W4 = 1 W2 = W3 = 1 + s. 



Equilibria with symmetric selection 

We start with only the modifier allele M present in the population so that Xi = x\ = 0 
for i = 5,6, 7, 8. Let the two-phase transformation of the frequencies, determined by (S2.3), 
be x' = T • Tx. This can be represented in two stages as 

x' = Tx, x = Tx, (52.21) 



where Xi = x; = x\ = 0 for i = 5, 6, 7, 8 and 

wxi = Wi(xi TnD) 
wx\ = Wi{xi Tr\D) 



1,2,3,4. (52.22) 



Here 



= W i X i ~ flD(wi -W2-W3 + W4) 

i=l 
4 

= ^2 - r x D{w x -W2-w z + w A ) (52.23) 



w 

4 

w 

D = X1X4 — X2X3 



D = X1X4 — X2X3. 
As wi = W4, W2 = W3, and w\ = W4, W2 = W3, we have 

, , _ WlWi 
WW 

, / W 2 W2, . 

x 2 ~ x 3 = —{X2 ~ X3). 

WW 
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(52.24) 



Hence 



wx = W\{x — r\D) = w 2 (l — ri)x + ~^W2- 



_ _ wi(l - r{)x + ^wi 

2(u>l - W 2 )(l ~ n)x + ^! + (1 - W 2 

w 2 (l — ri)x + 

x = 



(52.25) 



At equilibrium when x' = x, equations (S2.24) are satisfied by x\ = X4, X2 = X3, to- 
gether with x\ = X4, x 2 = xs- Hence, we look for symmetric equilibria of the form 
(#i, X2, xs, X4) = (x, I — x, \ — x, x). We obtain the following result: 

Result 2. On the boundary where only allele M is present, a unique symmetric equilib- 
rium x* exists where x* = (x*, | — x* , \ — x*,x*) and x* is the unique positive root of 
Q{x) = 0, where 0 < x* < \ and 

Q(x) = Ax 2 + Bx + C, 

A = n(l — ri)(wi + W2)(wi — W2) 

B = ^n(l - ri)(wi + w 2 )(wi - w 2 ) + ^ri (1 - (»i + u> 2 ) 2 
C = -^ri (1 - U> 2 (u>l + w 2 ). 

Proof of Result 2. 

If x = (x, I — x, I — x, x) , then x = (x, | — x, \ — x, x) and 

2 n V 1 ~ 2/1 "\ 2 1 

D = x — x =x , D = x — x =x 

V2 ; 4' V2 y 4 

w = 2(wi - w 2 )(l - ri)x + + (l _ y) w 2 

w = 2{w 2 - wi){\ - ri)x + —w 2 + [1 - —Jw 1 , 

since w\ = W2, W2 = w±. From (S2.22) we have 

wx = w\(x — r±D) = wi(l — ri)x + ~^w\ 



(52.26) 



(52.27) 



(52.28) 



2(w 2 - w x ){\ - ri)x + ^w 2 + (l-^)wi 
In (S2.28), substituting x in the expression for x, gives the quadratic equation Q(x) = 0 
given in (S2.25). Now as 0 < n < 1, 

Q(0) = (l - y) zo 2 (zoi + zo 2 ) > 0 (52.29) 
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and 



Q (I) = i 1 - 'i) + w i) > °- ( 52 - 3 °) 

Therefore Q(x) = 0 has a solution 0 < x* < |. It is the unique positive solution of 
Q(x) = 0 as Q (±oo) > 0 when w\ > w 2 and r > 0 so the other root of = 0 is 

negative. 

We next check the internal stability of x*. 

Result 3. The unique symmetric equilibrium x* is internally stable on the boundary 
where only the M modifier allele is present. 



Proof of Result 3. 

As x\ = X4 and X2 = x^ in x*, and in view of (S2.24), x* is internally stable if 
ww > wiWi and ww — W2W2, namely if ww > w\W2- But 



Therefore 



w = 2(wi - w 2 )(l - ri)x + W2 + —(wi - W2), 

W = 2(w 2 - Wi)(l - ri)X + Wi + "^"(^2 - Wi). 



ww = 2(w2 — W\){1 — ri)wx + w wi + —{w2 — w{) 



As wx = wi(x — r\D), we have 

ww — 2(w2 — Wi)(l — r\)w\X — 2r 1 (l — r 1 )(w2 — W\)w\D 



+ 



Using D = x 



2(wi - w 2 )(l - ri)x + w 2 + —(w 1 - w 2 ) wi + - wi) 



4' 



ww = ri(wi — W2) \ (1 — ri)x(w\ + w 2 ) + 



1 m 



Therefore as 0 < r± < 1 and W\ > W2, ww > w±w 2 if and only if 



(1 - n)x(wi + w 2 ) + 



(Wi + W 2 ) - W 2 



> 0. 



(52.31) 



(52.32) 



(52.33) 



(wi+w 2 )-W2 >+wiW2- (52.34) 



Treat the left side of the inequality (S2.35) as a function /() of ri, for 0 < r 1 < 
x,wi,W2 are parameters. In fact, 



/(I) 



1 



(Wi + W 2 ) - W2 



= |(u>i - w 2 ) > 0, 



/(0) = x(wi + w 2 ) - -w 2 . 



(52.35) 
1, where 

(52.36) 
(52.37) 
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So (S2.35) is satisfied, and x* is internally stable, if x* > 2 (w™+w 2 ) ' w i = 1 + s i 

W\ = 1. Thus we need x* > 2 ( s +2) ■ From (S2.28) we can write 



Q(x) = n(s + 2) |(1 -n)sx 2 
and we have Q(0) < 0 and 



i<i-n). + i(i-2)H.-i(i-S)} 



Q(^)=n( S + 2) -' 1 -;']y i ' <o. 

Since = 0, we must have x* > 2 ( s +2) as desired. 



(52.38) 



(52.39) 



Remark. Using (S2.38) it is easily seen that when s > 0 



Q(i)=ri(a + 2).->0. 



(52.40) 



Therefore as Q ^ 2 (s+2) ) < ^' we conc l u< ie that when w± > W2, 2 (s+2) < x * < V 



In 



addition, as D* = x* — \, at equilibrium D* < 0. 



External stability of the symmetric equilibrium 

We next check the external stability of x* to the introduction of the modifier allele 
m at the modifier locus. As x* is internally stable, its external stability is determined by 
the 4x4 matrix L* x , which is the linear approximation of the transformation x' = ToTx 
near x* restricted to the components x 5 ,xq,xj,xs of x. As our transformation is T o T, 
then L* x = L oL*, where we have 



L* = 



L = 





(a 


e 


h 


d\ 


1 


b 


/ 


g 


c 


w L 


c 


9 


f 


b 




\d 


h 


e 


a) 




/a 


e 


h 


d\ 


1 


b 


f 


9 


c 


W L 


c 


9 


f 


b 




u 


h 


e 


~a) 



(41) 



(42) 
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with 

w L = 2w 2 + (wi - w 2 ) [n + 4(1 - ri)x*\ 

(52.43) 

wl = 2zoi + (w 2 - wi) [ri + 4(1 - n)x] , 

and 

a = w 1 [2-r 2 - 2R(1 - r 2 )(l - x*)] 
b = W2 [R + r 2 - 2Rr 2 - 2R(1 - r 2 )x*} 
c = w 2 R[l - 2x*(l - r 2 )] 
d = 2w lX *R(l - r 2 ) 

(52.44) 

e = w 1 \(l-R)r 2 + 2R(l-r 2 )x*} 

f = W2 [ 2 -r 2 -R(l-r 2 )(l + 2x*)} 
g = w 2 R(l-r 2 )(l-2x*) 

h = Wl R[r 2 + 2(1 -r 2 )x*} 
a = w 2 [2-r 2 - 2R(1 - r 2 )(l - x)] 

b = W! [R + r 2 - 2Rr 2 - 2R(1 - r 2 )x] 

c = W!R[1 - 2x(l - r 2 )\ 

d = 2w 2 xR(l - r 2 ) 

(52.45) 

e = w 2 [(l - R)r 2 + 2R(1 - r 2 )x\ 
f = Wl [2 - r 2 - R(l - r 2 )(l + 2x)] 
~g = w 1 R(l-r 2 )(l-2x) 

h = w 2 R[r 2 + 2(1 - r 2 )x]. 

— * — * 

Here x = Tx* and x* = Tx. As L* x = L oL*, due to the special form of L* and L . 

L* x also has the same form, and the characteristic polynomial Q(X) of L* x factors into 

two quadratic polynomials, namely 

Q(\) = Q 1 (\)-Q 2 (\). (52.46) 

The analysis of Qi(X) and Q 2 (X) gives the following result. 

Result 4. The symmetric equilibrium x* is externally stable when n > r 2 and unstable 
when n < r 2 . Thus in this symmetric environment which changes every generation, larger 
recombination rates are favored, and the evolutionary stable recombination rate is 1. 
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Proof of Result 4. 

We outline here the steps of the proof. 

Let Q(A) be the fourth degree characteristic polynomial determining the external 
stability, with Q(X) = Q(\; r±, r 2 , R)- Then 

1. Q(A)=Qi(A)-Q 2 (A). 

2- Qi(A) = A 2 + Pi\ + qi for % = 1,2, where 

Pi > 0, 0 < qi < 1, Ai = Pi - 4qi > 0. 

3. Qi(X), for z = 1,2, has two positive roots whose product is less than 1. Both roots 
are less than 1 if and only if Qi(l) > 0. 

4. Qi(l) = Qi(l; r 1 ,r 2l R) = (ri-r 2 )/i(ri,r 2 , -R), where /i(ri,r 2 , -R) > 0 for all possible 
values of ri, r 2 , i?. 

5. When ri > r 2 , then Qi(l) > 0, and the two roots of Qi(\) = 0 are positive and less 
than 1. When r\ < r 2 , then Qi(l) < 0, and Qi(X) = 0 has a root larger than 1. 

6. When R = 0, then Qi(A;i?) = Q 2 (A;i?), and also Qi(l;i? = 0) = Q 2 (l;i? = 0). 

7. Q 2 (l; i?) is an increasing function of R. 

8. When R = 0, we have Qi(l; -R = 0) > 0 when n > r 2 , therefore Q 2 (l; R) > 0 for all 

when rx > r 2 . 

9. When ri > r 2 , both Qi(l) and (^2(1) are positive, hence all four roots of Q(\) = 0 
are positive and less than 1 making the equilibrium externally stable. 

10. When r\ < r 2 , then Qi(l) < 0, and so Q(A) = 0 has a root that is positive and larger 
than 1, and our equilibrium is externally unstable. 

All the steps above have been checked either analytically or using Mathematica for 
some computations, especially for those determining the sign of an expression, given the 
parameter space. Moreover, arguments involving the second polynomial Q2W are based 
on its following properties: 
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The second polynomial Q 2 (A) is the characteristic polynomial of 

1-9 b-c\ ff-g b-c\ 
e — h a — d J \ e — h a — d J 



T * — 

^ex — 



where 



w(a — d) = W\A w(a — d\ = W2A 
w(f-g)=w 2 A w(f-g)=w 1 A 
w(b — c) = W2B w(b — c) = w\B 
w(e — h) = W\B w(e — h) = W2B, 
A = 2 - r 2 - 2R(1 - n) = 2(1 - R) - r 2 (l - 2R) 

B = r 2 (l-2R). 



In fact, 



Thus 



1 / W1W2A 2 + w\B 2 [w\W2 + w 2 ) AB 
w lw l \ (w lW2 + w \) AB Wl w 2 A 2 + w\B 2 



w L w L (w L w L ) 
A 2 -B 2 = (A- B)(A + B) = 4(1 - R)[l - r 2 - R(l - 2r 2 )] 

g 2 (i) = 1 - 2wiW2A2 + ( " ? + ^ 32 + " ?W22 2 U 2 - £ 2 ) 2 . 

Let Q 2 (l) = Q 2 (l;i?). We compute using 

ggg(l) 2w lW2 ■ 2^-2(1 - r 2 )] + (w 2 + w|) ■ 2ff [-2r 2 ] 
<9i? wlwl 

+ • 2(A 2 - 5 2 ){2A[-2(1 - r 2 )] - 2S[-2r 2 ]}. 

(w L w L ) 



Hence, 



dQ 2 (l) _ 8w lW 2A(l - r 2 ) + 4 (u> 2 + w|) Br 2 
<9i? wlwl 

-^%(A 2 -B 2 )[A(l-r 2 )-Br2], 
[w L w L ) 

^Mll > 0 « (w L w L ) [8 Wl w 2 A(l - r 2 ) + 4 (u> 2 + u> 2 ) Sr 2 ] 

> 8zo 2 «;i(A 2 - 5 2 )L4(1 - r 2 ) - Br 2 ], 
A 2 - £? 2 = 4(1 - r)[l - r 2 - i?(l - 2r 2 )]. 
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Thus > 0 if 

w L w L A(l - r 2 ) > 4wiiy 2 [l - r 2 - R(l - 2r 2 )]A(l - r 2 ). 
As A > 0 and 0 < r 2 < 1, this is true if 

wlwl > 4u>iu> 2 [l - r 2 - i?(l - 2r 2 )]. 

But wlWl > ^wiw 2 and 1 — r 2 — i?(l — 2r 2 ) < 1, and the above holds. 

This sums up to the fact that Q 2 (1;R) is monotone increasing in R, thus proving 
point 7. But since Q 2 (l, 0) > 0 when r\ > r 2 , therefore Q 2 (l; R) > 0 when r\ > r 2 for all 
0 < R < 1. 
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